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Abstract 

The analytic solutions of the one-dimensional Schrodinger equation for the trigonometric Rosen- 
Morse potential reported in the literature rely upon the Jacobi polynomials with complex indices 
and complex arguments. We first draw attention to the fact that the complex Jacobi polynomials 
have non-trivial orthogonality properties which make them uncomfortable for physics applications. 
Instead we here solve above equation in terms of real orthogonal polynomials. The new solutions 
are used in the construction of the quantum-mechanical superpotential. 



PACS numbers: 02.30.Gp, 03.65.Ge, 12.60.Jv 



I. INTRODUCTION 

n 

Supersymmetric quantum mechanics was originally proposed by Witten |l| as a simple 
learning ground for the basic concepts of supersymmetric quantum field theories but soon 
after it evolved to a research field on its own rights. Supersymmetric quantum mechanics 
starts with the factorization of one-dimensional Hamiltonians, 

H ^ = ~^ + v ^ « 

according to H(y) = A + (y)A~(y) + e with A ± (y) = (^±-^=-^ + U(y) s j where U(y) is the 
superpotential. Next it proves that if ip n (y) is an exact solution of the H(y) eigenvalue prob- 
lem, H(y)i/; n (y) = E n i[) n (y), then A(y)ij; n (y) is an eigenfunction to H(y) = A~{y)A + (y) + e 
corresponding to same eigenvalue. In other words, knowing the superpotential allows to 
generate the H (y) spectrum from the spectrum of H (y) and vise versa. Moreover, in case of 
zero ground state (gst) energy, knowing U(y) allows to solve A~ (y)ipg S ^(y) = and obtain 
the ground state wave function. In other words, knowing the ground state wave function 
allows to recover the superpotential as 

At that stage one uses the isospectral pair H(y), and H(y) in the construction of the super- 
Hamiltonian H(y) as H(y)=diag( H(y), H(y) ) and upon introducing "charges" as 




Q(y) = I A ° {y) ° I , Q'(y)= I I *" I , (3) 

proves them to satisfy the following algebra: 

{Q(y), Q\y)} = H(y), {Q(y), Q(y)} = {Q\y), Q\y)} = 0, 

[Q(y),H(y)] = [g t (y),H(y)] = 0. (4) 

The relationship to the field-theoretic SUSY is then established through the observation 
that in case A~(y)ipg S ^(y) ^ then the charges do not annihilate the respecive vacua, 
Q(y)i/jgst(y) 7^ 0, and Q^iy^gstiy) 0> an< i SUSY is spontaneously broken. On the 
contrary, when A~ (y)ipg S ^(y) = 0, then the charges annihilate the ground states which is 
equivalent to the absence of charge condensates there, and thereby to SUSY realization 



in the manifest (multiplet) Wigner-Weyl mode. In this manner supersymmetric quantum 
mechanics relates to SUSY in field theory where the role of Q(y) and Q'(y) is taken by 
boson-fermion (and vise versa) ladder operators. 

Above considerations clearly reveal importance of knowing the exact solutions of the quan- 
tum mechanics Hamiltonians. These solutions are furthermore important in the construction 
of higher dimensional charge algebras with more but two Hamiltonians (hierarchy of Hamil- 
tonians) |2|. The supersymmetric quantum mechanics manages a family of exactly soluble 
potentials |3j, one of them being the trigonometric Rosen- Morse potential (tRMP). As long 



as this potential is obtained from the Eckart potential 4] by complexification of the argument 
and one of the constants, also its solutions have been concluded from those of the Eckart po- 
tential by same procedure. In so doing, one ends up with Jacobi polynomials with complex 
indices and complex arguments. However, such complex polynomials are not comfortable 
for physical applications basically in view of their non-trivial orthogonality properties 0, f| . 

We here make the case that the trigonometric Rosen-Morse potential is exactly 
soluble in terms of a family of real orthogonal polynomials and present the 
solutions. 

The paper is organized as follows. In the next Section we present the tRMP derivation 
form the Eckart potential and draw attention to the non-trivial orthogonality properties of 
the Jacobi polynomials with complex parameters and arguments. In Section III we solve 
analytically the one-dimensional Schrodinger equation with the trigonometric Rosen-Morse 
potential and present the solutions. In Section IV we employ the exact ground state wave 
function in the construction of the tRMP superpotential. The paper closes with brief con- 
cluding remarks. 



II. THE TRIGONOMETRIC ROSEN-MORSE POTENTIAL AS COMPLEXIFIED 
ECKART POTENTIAL. 

Before proceeding further we first introduce a properly chosen length scale d and change 
variables in the one-dimensional Schrodinger equation (JIJ) to dimensionless ones according 
to 

z=2> v(z) = V(dz)/(h 2 /2md 2 ), e n = E n /(h 2 /2md 2 ) . (5) 
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Next we employ the Eckart potential 

v(z) = — 26coth2; + a(a — 1) csch 2 z , 
where b > a 2 . The exact solutions to the Eckart potential read: 



(6) 



b 

x = coth Z , Pn = — — • (7) 



n + a 



Here p^j 3n ~ n ~ a -(Pn+n+a)) 



(x) are the well known Jacobi polynomials with n < (b 1 ^ 2 



and Cfi IS ct normalization constant. Equation (JZj) equivalently rewrites to 

ij) n (x) = c n (x 2 — l)~( n+a )/ 2 e _/3n arccoth x p(/ 3 ™- n - a >-(/ 3 ™+ n + a ))( x ) 

The corresponding energy value spectrum is determined by 

^2 b 2 
t n = -(n + a - 



(8) 



(9) 



(n + a) 2 

Form and energy spectrum of the Eckart potential are illustrated by Fig. I. Let's now 
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Fig. I. Eckart potential. The solid curve represents the potential while the dashed lines are the 
energy levels. Notice that the argument of this potential is unbound from above, i.e. < z < oo 
and the number of bound states limited. The parameters of the displayed potential take the 

values a = —1, b = 50. 

complexify the argument of the Eckart potential and one of its constants according to 

z — > — i z , or, equivalently, x — > ix; b — > ib. (10) 
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Substitution of Eq. (jTUJl into (JEJ) results in 

v (z) = —2b cot z + a(a — 1) esc 2 z , (11) 

and thereby in the trigonometric Rosen- Morse potential [3] shown on Fig. [H] 
In the literature jsj the solution of the Schrodinger equation with the trigonometric Rosen- 



v(z) 




Fig. II. The trigonometric Rosen-Morse potential. The solid line represents the potential, while 
the dashed lines are the energy levels. Notice that the argument of this potential is bound 
between < z < n and the number of states is unlimited. The potential parameters take the 

values a = 1, b = 50. 

Morse potential is concluded from Eq. (jHJ) through complexification of b and x leading to 

ll> n {ix) = C n ((ix) 2 - l)-("+ a )/ 2 e- ftl arccoth ix Pf Pn-n-a)-(if3 n +n+a))( ix j _ ^ 

In other words, these solutions need the Jacobi polynomials with complex indices and com- 
plex arguments. Unfortunately, the complex Jacobi polynomials are not very well suited for 
physics applications. Suffices to write down these polynomials explicitly, 

«*) = ( JT + + B %„ (-•" + A + B + 1; A + 

( „ ,l-ix\ T(A + 1) 

2 F 1 [-n,n + A + B + l;A + l\—- ' 



T(n + A + B + l)r(-(n + B)) 
J 1 t n+A+B {l - t)-(-+ B + 1 ) (\ - J-^j n dt , 
A = i/3 n — n — a , B — —(i(3 n + n + a) , (13) 



where 2-^1(0, b; c\x) is the well known hypergeometric function, and (...)„ is the Pochhammer 
symbol of the expression in the parentheses, in order to become aware of the calculational 
difficulties to be expected. It is obvious that one has to worry about the interplay between 
indices and integration contours, a subject studied in 0, |(| . There the authors claim de- 
pendence of the orthogonality properties on the indices and the integration contours. In 
order to avoid all those difficulties we here search for real solutions of the one-dimensional 
Schrodinger equation with the trigonometric Rosen-Morse potential. 

III. THE ONE-DIMENSIONAL SCHRODINGER EQUATION FOR THE 
TRIGONOMETRIC ROSEN-MORSE POTENTIAL. 

A. The Sturm-Liouville method and the Rodrigues formula. 

For the sake of self-sufficiency of the representation we here review in brief the basics of 
the Sturm-Liuoville technique for solving second order differential equations. 
The method of Sturm-Liouville applies to differential equations of the type 

-j— (pWt^I + q{x) y = Xw(x) y , (14) 
a x \ a x J 

and searches for the A values that allow for a solution. The solution of Eq. (fT4*j) . where 
w(x) stands for the "weight", or, "density" function, are eigenf unctions of a Hermitian 
differential operator on the space of functions defined by the initially conditions. Some 
special cases of the Sturm-Liouville type of differential equations allow for solutions by 
means of the Rodrigues formula pj. To be specific, for q(x) = 0, and p(x) = w(x)s(x), 
where s(x) at most a second order polynomial, the solution of Eq. ()14|) is given by a family 
of orthogonal polynomials. The classical polynomials of Hermite, Laguerre, Legendre and 
Jacobi are prominent examples for that. 



In order to create the orthogonal polynomial solutions one considers the function 
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1 d m 

C m (x) = J w(x) s(xD , (15) 

w{x) a x m 



with Ci(x), w(x), and s(x) satisfying the following conditions: 
1. C\(x) is a polynomial of first order, 



2. s(x) is a polynomial of at most second order and real roots, 

3. w(x) is real, positive and integrable within a given interval [a, b], and satisfies the 
boundary conditions 



w(a)s(a) = w(b)s(b) = . 



(16) 



Above three conditions seem quite restrictive indeed but nonetheless they allow for the 
construction of all the classical orthogonal polynomials reported in the standard textbooks 
0, 0], 0. The table shows all the ingredients of the Rodrigues formula required for the 
construction of the respective orthogonal polynomials. 



TABLE I: Special functions and their characteristics. 



Name 


Symbol 


w{x) 


s(x) 


Interval 


Conditions 


Hermite 


H m (x) 


e -x 2 


1 


( — OO, Cxd) 




Laguerre 




x v e~ x 


X 


[0,oo) 


(u>-l) 


Jacobi 




{l-x) u {l + xy 


(1-x 2 ) 


[-1,1] 


(v,H> -1) 


Gegenbauer 


C rn (x) 


{1 _ x 2 ) X-l/2 


(1-x 2 ) 


[-1,1] 


(A > -1/2) 


Legendre 


P m {x) 


1 


(1-x 2 ) 


[-1,1] 




Chebyshev, Type I 


T m (x) 


(1-X 2 )-V2 


(1-x 2 ) 


[-1,1] 




Chebyshev, Type II 


U m (x) 


(1-x 2 ) 1 / 2 


(1-x 2 ) 


[-1,1] 




This work 




+ e -^ cot " la; 


(1 + x 2 ) 


( — CXD, Cxd) 


(pL = m + 1 + a) 



If in addition one demands an orthonormalized set of polynomials, one has to introduce 
into Eq. (|15|) an additional constant, here denoted by K m , according to 

1 d ri 



C m (x) 



K m w(x) d x r 



[w(x) s(x) m ) 



In terms of C m (x), the Sturm-Liouville equation takes the form 

(w(x) s(x) ^ G ™ X ) = -X m w(x) C m (x) 
a x \ ax) 

or, equivalent ly, 



s[x) 



d 2 C m (x) 
d x 2 



1 



W(X) 



d s(x)w(x)\ d C m (x) 



d x 



d x 



\- X m C m (x 



(17) 



(18) 



(19) 



/ , r d CAx) 1 . „ .d 2 s(x) \ 



where 



Notice that 



K 1 w(x) \ d x 

All classical polynomials can be obtained through above procedure and vise versa. The 
polynomials that can be obtained from the Rodrigues formula and which satisfy the three 
conditions mentioned above are necessarily the classical orthogonal polynomials, a result 
due to Ref. Q. 



CxV) = T-Vt ( ~ S{X J WiX) ) • (21) 



B. Solving Schrodinger's equation for the trigonometric Rosen-Morse potential. 

In this Section we present the solution of the one-dimensional Schrodinger equation for 
the trigonometric Rosen-Morse potential as obtained in Ref. ^Jj and without any reliance 
on the complex Jacobi polynomials. For this purpose we first have to reshape Schrodinger's 
equation 

d 2 R(z) 

— — ^ + (2b cot z-a(a + 1) esc 2 z + e) R(z) = , (22) 
a z l 

to the Sturm-Liouville form in Eq. (|19p. To do so we first change variables to 

R( z ) = e - az/2 F{z) , (23) 



with a being a constant and then substitute in Eq. (J22J). After some simple algebraic 
manipulations one arrives at 

d z 2 

Changing once more variables to x — cotz in which case F(z) becomes a function of x 
denoted by f(x), i.e. 



a ^Tz^ + ( 2bcotz - a ( a + 1 ) csc2z + ((I)' + e )) F ( z ) = ° ■ ( 24 ) 



F(z) -> f(x) , x = cot z, (25) 

Eq. (|24)l takes the form 



^2bx - a{a + 1)(1 + x 2 ) + ((f )' + e ) 



f(x) + (1 + x 2 ) 



2\2^ f( X ) 



d x 2 

d f(x) 



+ 2 (i^)(f + .)^ = o. (a6) 



Finally, upon substituting the factorization ansatz 

f(x) = (l+x 2 )- {1 -V /2 C(x) , (27) 

into Eq. pfij) and a subsequent division by (1 + x 2 )( 1+ ^)/ 2 one finds as intermediate result 
the equation 



-a I 



l-P) + 2b)x+ ((f) 2 -(l-6) 2 + e) 

o. d 2 C(x) /a „ \ d C(x) 



(-/3(1 - /3) - a(a + 1)) + ^ | C(x 



(28) 

If Eq. (J2*B*|) is to coincide in form with Eq. ([19)1. following conditions have to be fulfilled: 

-a(l-/3) + 26 = 0, (29) 
(|) 2 -(l-/3) 2 + e = 0. (30) 

Substitution of Eqs. flU, fl3H into Eq. J2SD amounts to 

(1 + x2) ^M + 2 (| + ^ d_C^) + _ p) _ ^ + 1))c{x) = Q (31) 

Equation (jSU relates to Eqs. (fTTIil-CT via 

- - /?) - a(a + 1) = -m(2/3 + m - 1) . (32) 

In now determining /3 from Eq. (jS2J), substituting it into Eqs. (|2~H|) . (jSHJ), and shifting m to 
to — > n — 1, the following n dependent constants are found: 

2b 

Pn = -{n + a) + 1 , a n = — ■ — , (33) 

n + a 

b 2 

e n = (n + a) 2 - , (34) 

(n + a) z 

with n > 1. 

In this way one encounters ry(x) and s(x) as 

Wn (x) = (1 + x 2 )- (n+a) e- Q " arccot x , (35) 
s{x) = 1+x 2 . (36) 
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The polynomials which resolve the tRMP Schrodinger equation are now obtained in exploit- 
ing the Rodrigues formula (fTTjl when rewritten in terms of n as 

1 d^ 1 

The lowest Cn ,b \x) polynomials obtained in this fashion read 
C[ a ' b \x) = ±- , (38) 



(39) 



Ct b \x) = J- ((1 + a)(2a + 3)x 2 - 2(2a + 3)JL-x + (t^j " ^ + a ) ) ) ' ( 40 ) 

Cf' fe) (x) = -jj- ( -(! + a )( 2a + 3) (2 + a)x 3 + 3(a + 2) (2a + 3) ^ . x 2 
A 4 \ (4 + a) 

-3(2 + a) ( V & - (1 + a)] x + ( ( - (3a + 4)—^— 1 ) . ( II ) 



(4 + a) 2 v 7 V( 4 + a ) 3 4 + a 

Ct' b \x) = 4r ((l + a)(2a + 3)(2 + a)(2a + 5)x 4 -4(2a + 3)(2 + a)(2a + 5) - 6 . x 3 
K 5 \ (5 + a J 

+6(2 + a) (2a + 5) - (1 + a) ) x 2 

V(5 + a) 2 J 

(2b 3 b \ 

- 4 ( 2o + 5 '((5T^-( 3o + 4 W a: 

/ 46 4 46 2 \\ 

+ [WTar " (FT^ (3 ° + 5) + 3(2 + '( 1 + 'j) ' < 42 > 

where x = cot z. 

Above polynomials solve exactly Eq. (J31|) which can be immediately cross-checked by back- 
substituting Eqs. ()42j) into Eq. (|31|). Employing symbolic mathematical programs is quite 
useful in that regard. 

Notice that the solution was found under less rigid requirements but the ones 
listed immediately after Eq. ([T5]) above. Indeed, 

• the roots of our s(x) function are imaginary (this is the only place where 
the complexification of the Eckart potential seems to have left footprints), 

• Equation (|TH|) turned out to be a rule that allows to be broken. 

In that sense, the orthogonal functions found here must belong to a different 
class of orthogonal functions. 
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Their orthogonality is obtained as 

"°° dx 



(w n (x) ) 1/2 C^' b) (w n , (x) ) 1/2 C { *' b) = 5 n n , . (43) 

S{X) 

Equation (}4*3*|) shows convincingly that the new solutions have well defined orthogonality 
properties on the real axes, which qualifies them as comfortable wave functions in quantum 
mechanics applications. 

The orthogonality condition for the wave functions R n (z) reads 

/ dz R n (z)(R n ,(z)y = 5 nn ,. (44) 
J o 

A version interesting for physical application (see concluding Section) is the one with a = 0, 

v(z) = -26 cot z, (45) 
in which case the normalization constant is calculated as 

^( {n[)2 :l- e r /n) V /2 - (46) 



(47) 



Ab{b 2 + n 4 ) 

The associated energy spectrum is given then by 

b 2 



e n = n 2 



Correspondingly, the wave functions for this case simplify and are found as 

R x (z) = e~ bz smz C{°' b \cotz) , 
R 2 (z) = e~ bz/2 sin 2 zCi°' b) (cot z) , 

R n (z) = e- 62/n sin"z C(°' 6) (cotz) . (48) 
Wave functions for the first two (unnormalized) levels with a ^ are displayed in Figs. III. 



IV. THE TRIGONOMETRIC ROSEN-MORSE SUPERPOTENTIAL. 

In this Section we derive the trigonometric Rosen-Morse superpotential from the exact 
ground state solution which for a ^ has been calculated as 

Ri(z) oc e- bz/{a+1) sm a+1 z . (49) 

11 





Fig. III. Wave functions for the first two levels in the trigonometric Rosen-Morse potential for 

a =0.25 and b =1. 

Substitution of the latter equation (jlHJ) into Eq. (J2J) amounts to the following superpotential, 

b 



U{z) 



a + 1 



+ (a + 1) cot z 



shown in Fig II VI Correspondingly, the A (z) operators are obtained as 



(50) 




Fig. IV. The trigonometric Rosen-Morse superpotential. The values of the displayed 
superpotential parameters are same as in Fig. [Hj 

A ± (z) = ±— + (a+ l)cotz — . (51) 

w dz a + 1 v ' 

The corresponding Hamiltonian is then obtained identical to Eq. as should be. The 

super symmetric companion of H(z), which is H(z), becomes 

d 2 

H(z) = -— - 26 cot 2 + (a + l)(a + 2) esc 2 ^ , (52) 
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and one finds the supersymmetric companions to the solutions as 

Rn(z)\(a,b) = Rn-l(z)\(a+l,b) , (53) 

where n > 1, and R n ~i(z)\( a +i,b) are among the exact solutions of the trigonometric Rosen- 
Morse potential. 

As a further possible application of the solutions found here we wish to mention the con- 
struction of the so called hierarchy of Hamiltonians |2j, where one needs to have at ones 
disposal exact orthonormalized functions for all levels because in this case one can pick up 
any energy level, e n , and its wave function R n (z). Obviously Eq. (|4*Kj) fully qualify for that 
purpose. 



V. CONCLUDING REMARKS. 



To recapitulate, the knowledge on the real exact solutions of the one-dimensional 
Schrodinger equation for the trigonometric Rosen-Morse potential obtained in this work 
allows for a straightforward construction of all the necessary ingredients of the supersym- 
metric quantum mechanics. Compared to the text-book complex solutions, the exact real 
solutions reported here have the advantage to significantly simplify the calculations. 

In addition, the trigonometric Rosen-Morse potential is more but just one of the few 
exactly soluble quantum mechanical potentials. It joins as a new member the smaller and 
important group of potentials which generate physically relevant spectra. Indeed, in particle 
physics one encounters [l2j the measured excitation spectra of the nucleon (A), and the so 
called A particle, in turn displayed in Fig. |V| 

According to an observation due to Refs. [lflj . [l^ and references therein, those spectra 
repeat with an amazing accuracy the degeneracy patterns of the levels of the electron with 
spin in the hydrogen atom but are characterized by very different mass splittings. The 
following mass formula describes pretty well the averaged positions of the three narrow 
mass bands containing the series of (n — 1) parity doubled states (with n — 2,4, and 6) 
and with spins ranging from ^ to (n — |) together with the one unpaired maximal spin 
( n ~~ l) P s ^ a ^ e °f either natural, or, unnatural parity, 

n 2 — 1 1 

M (n;/) - M ( ° n;7) = gi — — - fj— , I = N,A. (54) 
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Fig. V. Experimentally observed baryon resonances (l.h.s.) iV and (r.h.s.) A. The dash-point 
lines represent the mass average. Notice that the resonances with masses above 2000 MeV are of 
significantly lower confidence but those with masses below 2000 MeV where the degeneracy is 
very well pronounced. Empty squares denote predicted ("missing") states. 

Comparison of the baryon mass formula with Eq. ()47j) reveals coincidence between the 
baryonic and the tRMP spectra. The underlying constituent dynamics of the excited baryons 
has been found to be that of a quark-di- quark system Q|. In fact, in order to describe the 
particle spectra one needs to solve the three dimensional Schrodinger equation but this does 
not cause much problems because upon separation of the variables in polar coordinates the 
radial part of the three dimensional Schrodinger equation reduces to that very same one- 
dimensional Schrodinger equation (|22|) up to the centrifugal barrier term, 1(1 + l)/z 2 (work 
in progress). All in all, the trigonometric Rosen-Morse potential and its real orthogonal 
polynomial solutions open new venues in the calculation of interesting observables in both 
super symmetric quantum mechanics and particle spectroscopy. 

VI. ACKNOWLEDGMENTS. 



We thank Dr. Alvaro Perez Raposo for his spontaneous and vivid interest in our results. 
Work supported by Consejo Nacional de Ciencia y Technologia (CONACyT) Mexico under 
grant number C01-39280. 



14 



[1] E. Witten, Nucl. Phys. B 188, 513 (1981). 

[2] C. V. Sukumar, J. Phys. A:Math. Gen. 18, 2917 (1998); 

C. V. Sukumar, AIP proceedings 744, eds. R. Bijker et al., Supersymmetries in physics and 
applications, (New York, 2005), p. 167. 

[3] F. Cooper, A. Khare, U. P. Sukhatme, Supersymmetry in Quantum Mechanics, (World Sci- 
entific, Singapore, 2001). 

[4] R. De, R. Dutt, U. Sukhatme, J. Phys. A: Math. Gen. 25, L843 (1992). 

[5] A. B. J. Kuijlaars, A. Martinez-Finkelshtein, R. Orive, Orthogonality of Jacobi Polynomials 
with General Parameters, E-Print ArXiv: math.CA/0301037 (2003). 

[6] B. Beckermann, J. Coussement, W. Van Asshe, E-Print ArXiv: |math.CA/0311055| (2003). 

[7] Phylippe Dennery, Andre Krzywicki, Mathematics for Physicists, (Dover, New York, 1996). 

[8] M. Abramowitz, I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs 
and Mathematical Tables, (Dover, 2nd edition, New York, 1972). 

[9] http://en.wikipedia.org/wiki/Main_Page 



[10] C. W. Cryer, Boll. Un. Mat. Ital. 25 (1), 1 (1970); 

Renato Alvarez Nodarse, Las funciones especiales y la fisica matematica (in Spanish), 
http:/ /euler. us.es/ renato/clases.html 

[11] Cliffor Benjamin Compean Jasso, Baryon spectra in a quark- di- quark model with the trigono- 
metric Rosen-Morse potential, (in Spanish), MS thesis, Institute of Physics, Autonomous 
University of San Luis Potosf, Mexico, September (2005). 

[12] L. Alvarez-Gaume et al, Phys. Lett. B, 592, 1 (2004). 

[13] M. Kirchbach, M. Moshinsky, Yu. F. Smirnov, Phys. Rev. D 64, 114005 (2001). 
[14] M. Kirchbach, Rev. Mex. Fis. 50, Supl. 2, 54 (2004). 



15 



